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ABSTRACT 

We  consider  general  linear  parabolic  equations  in  a given  time  dependent  domain 
and  we  describe  a general  class  of  Galerkin-type  approximations  which  are  continuous 
with  respect  to  the  space  variables,  but  which  admit  discontinuities  with  respect  to 
time  at  each  time  step.  Unconditional  stability  is  proved  and  a general  error  estimate 
is  established.  These  results  are  applied  to  certain  finite  element  methods  based  on 
space-time  finite  elements. 

AMS(MOS)  Subject  Classification  - 65N05 

Key  Words  - Galerkin-type  method,  Finite  elements,  Parabolic  equations,  Moving  boundary, 
Stability,  Error  estimate. 

Work  Unit  Number  7 - Numerical  analysis 

EXPLANATION 

Most  numerical  methods  for  parabolic  equations  are  based  on  a space  discretization 
which  is  independent  of  time  and  are  not  appropriate  for  problems  in  a variable  domain. 
These  problems  can  be  approximated  by  using  finite  elements  which  are  relative  to  both  the 
space  and  time  variables;  such  elements  have  been  used  for  free  boundary  problems  in  heat 
conduction  (Stefan  problem)  and  fluid  flow,  but  no  mathematical  results  have  been  proved. 
In  this  report,  we  consider  the  case  of  a given  moving  boundary  and  we  prove  stability 
and  convergence  with  error  estimates  for  a general  class  of  Galerkin-type  methods  which 
includes  the  case  of  space-time  finite  elements. 
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1.  Introduction 

Most  numerical  methods  for  parabolic  partial  differential  equations  are  based  on  a 
space  discretization  which  is  independent  of  time;  this  is  the  case  for  most  finite  ele- 
ment or  Galerkin-type  methods  which  have  been  studied  in  recent  years  ([2],  [6],  [7],  [21], 
[23],  [24],  [26],  [27]).  These  methods  are  not  appropriate  for  parabolic  problems  in  time- 
dependent  domains,  in  particular  for  time-dependent  free  boundary  problems. 

In  order  to  solve  such  problems,  numerical  methods  based  on  space-time  finite  ele- 
ments, i.e.  on  finite  elements  which  are  relative  to  both  the  space  and  time  variables 

have  been  proposed  and  numerically  tested  ([3],  [4],  [11],  [12]).  Let  us  notice  that 

(**) 

space-time  finite  elements  have  been  first  considered  by  J.  T.  Oden  [19]  , but  only  in 

the  case  when  the  elements  are  the  Cartesian  product  of  a space  element  by  a time  interval, 
which  yields  a space  discretization  which  remains  fixed  in  time. 

Another  finite  element  method  to  deal  with  variable  domain  has  been  studied  by 
M.  Mori  [16] , [17]  for  the  Stefan  problem;  it  is  based  on  finite  elements  in  space  which 
depend  continuously  on  time.  This  method  is  indeed  a particular  case  of  the  generalized 
Galerkin  method  studied  by  Mignot  [15];  in  the  same  paper,  Mignot  studies  other  methods 
for  parabolic  equations  in  a variable  domain:  method  of  fixed  auxiliary  domain,  method 
of  elliptic  regularization.  Let  us  finally  mention  the  finite  difference  methods  studied 
by  the  author  [9]  for  parabolic  equations  of  order  2 which  can  degenerate  or  admit  singu- 
larities at  the  initial  time  as  well  as  on  the  boundary  of  the  variable  domain. 

In  the  present  paper,  we  consider  general  parabolic  equations  in  a given  time- 
dependent  domain  and  we  present  a general  class  of  Galerkin-type  methods  with  applications 
to  space-time  finite  elements.  These  methods  are  different  from  those  which  have  been 

* Visiting  from  Centre  d'Etudes  de  Limeil,  Commissariat  & l'Energie  Atomique,  B.  P.  27, 
94190-Villeneuve  St.  Georges,  France. 
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studied  previously  (including  the  space-tune  finite  element  methods  of  [ 3 ] # (4] , [11], 

112]):  the  approximations  are  continuous  with  respect  to  the  space  variables  for  each 


fixed  time,  but  they  admit  discontinui t ies  with  respect  to  the  time  variable  at  each  time 
step;  in  particular,  the  elements  can  be  chosen  arbitrarily  at  each  time  step  with  no 
connection  with  the  elements  corresponding  to  the  previous  step. 

For  more  simplicity,  we  first  consider  a model  problem:  the  Dirichlet  problem  for 
the  heat  equation.  This  problem  is  presented  in  Section  2;  we  derive  an  integral  relation 
which  is  satisfied  by  the  solution  and  which  is  the  basis  of  the  numerical  method.  In 
Section  3,  we  describe  a general  class  of  time-discontinuous  Galerkin-type  approximations 
and  prove  unconditional  stability;  the  discontinuities  are  treated  as  in  the  discontinuous 
finite-element  methods  cf  Reed  and  Hill  [22]  and  Lesaint  and  Raviart  [13]  for  the  stationary 
neutron  transport  equation;  the  same  technique  has  also  been  used  by  Fini  [20],  In 
Section  4,  we  establish  a general  error  estimate.  In  Section  b,  we  apply  the  previous  re- 
sults to  space-time  finite  element  methods;  the  order  of  accuracy  can  bo  made  arbitrarily 
high  by  choosing  finite  elements  of  corresponding  order.  In  Section  6,  we  give  a simple 
example  in  the  one-dimensional  case;  we  give  the  explicit  form  of  the  discrete  equations 
for  rectangular  space-time  finite  elements  of  order  1,  which  yields  a finite  difference 
analogue  of  our  method.  Finally,  in  Section  7,  we  extend  the  method  and  the  previous  re- 
sults to  general  parabolic  equations  of  order  > 2 with  general  boundary  conditions. 


r 


2.  A model  problem 


Let  us  consider  a time  interval  0 t <_  T and  let  (l(t)  be  a bounded  domain  in 
lRm,  m positive  integer,  which  depends  continuously  on  t e [0,T]  Let  r (t)  be  the 
boundary  of  (5(t)  and 

A T = { (x , t ) ; x e fi( t)  , 0 < t < T } , 

l = { (x,t)  ; x e r (t)  , 0 < t < T } . 

is  a (m+1) -dimensional  domain  in  IR m * n + and  £T  is  its  "lateral"  boundary.  We 
assume  that  £ is  piecewise  smooth.  For  simplicity,  we  will  use  the  same  notation  ft(t) 
for  the  domain  (2(t)  c and  for  the  corresponding  section  of  A^,  i.e.  the  set 

{ (x , t)  ; x £ fi(t)  , t fixed}  . 

Let  f and  u'1  be  two  given  functions,  f c l?  (A  ) , uG  c L2(0(O)),  and  let  A 
denote  the  Laplacian  operator  with  respect  to  the  space  variables.  We  consider  the  problem 


a)  — - Au  = f in  A , 
1 


b)  u = 0 


c)  u = u 


on  ) , 

tT 


in  (1(0)  . 


For  the  existence  and  uniqueness  of  a classical  solution  u,  under  suitable  hypotheses 
on  f and  uG,  see  (81 , (9) ; for  the  existence  and  uniqueness  of  a weak  solution,  see 
(14),  [15).  In  this  paper,  we  will  assume  that  the  solution  u is  sufficiently  regular 
for  the  validity  of  the  error  estimates. 

We  will  use  the  following  notation: 
d = Max{diam  (l(t)  ; 0 < t < T]  , 


<",")(l(t)  = inner  product  in  L (fl(t))  , 

I * Ijj  (t)  = n0rm  in  b2((l(t))  , 

((•»•))_  = inner  product  in  1.2(G),  where  G is  a subdomain  of  A , 

G T 

II  • ||  „ = norm  in  I,2(G)  . 

G 

The  same  notation  will  be  used  for  vector  valued  functions  in  (L2(G))IT1  . Thus, 


((grad  if/,  grad  ^))  = //  grad  \);  * grad  dxdt  , 

G G 
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where  i>  and  f axe  two  arbitrary  smooth  enough  functions. 


Now,  let  and  be  two  arbitrary  numbers  such  that  0 < tq  < Ti  <_  T and  let 

G « G(tq,t^)  denote  the  intersection  of  J>  ^ with  the  strip  < t < t^,  i.e. 

G « {(x,t);  x £ ait) , < t < t^}  . Consider  the  bilinear  form 

(2.2)  8 = -(  ^'ft'  )G  + ((9rad  9tad  <f))Q  + j 

Let  #(G)  be  the  space  of  all  Lipschitr-continuous  functions  v>  defined  on  G (closure  of 
G)  and  which  vanish  on  the  lateral  boundary  of  G,  i.e.  on  I n G . Then,  a classical 
integration  by  parts  shows  that  the  solution  u of  problem  (2.1)  satisfies  the  integral 

relation 

(2.3)  £ (u,*>)  - ((f,*)>  for  all  <p  £ *(G)  and  for  all  G = G(t  ,r  ) 

G G U i 

vith  0 <_  tq  < t <_  T 

This  relation  is  the  basis  of  the  numerical  method  described  in  the  following  section. 


-f*)  It  is  possible  to  relax  the  hypotheses  on  the  function  f in  relation  with  the 

regularity  of  u . This  question  is  considered  in  Section  fc  where  a precise  vari- 
ational formulation  is  e*  n for  more  general  parabolic  problems. 


3.  Discontinuous  approximations 


Let  {t  ; 0 <_  n <_  N)  be  a finite  sequence  of  real  numbers  with  t°  = 0,  t°  < tn+* 


and  t - T . Let  fj"  - (5(tn),  g"  - G(tn,t"+1)  and  g"  - G n - fi n 
tn  < t t"  (see  Figure  1). 


{ (x,t) ; x £ fl(t)  , 


Figure  1 : The  discretization  of  the  domain  with 

respect  to  time  (one-dimensional  case) 


Let  be  a finite  dimensional  subspace  of  4>(Gn),  for  0 < n < N-l.  and  let  V.  be 

h — — h 

the  space  of  all  functions  v.  defined  n J>  such  that  their  restriction  to  each  G° 

h T 

coincides  with  the  restriction  to  Gn  of  a function  £ $"  . The  functions  e V. 

h h h h 

are  in  general  discontinuous  at  the  time  t = tn;  we  will  use  the  notation  v”  = v^("(tn) 

h h 

for  0 < n < N and  = lim{v  (.,t n+e)  ; c > 0,  e ■+■  0}  for  0 < n < N-l;  it  follows 

— — h h — ~ 

from  the  definition  of  the  space  V.  that  we  have  v"  = lim{ v.  ( * ,tn-t)  ; z > 0,  z 0} 

h h h 

for  1 £ n 1.  N • 

We  approximate  problem  (2.1)  by  the  following  problem  which  is  a discrete  analogue 

of  (2.3).  Discrete  problem:  Find  u.  £ V.  such  that  u°  = u°  and 
n n n 

(3-1)  * nW  " <(f'V>n  ’ 

U (r 


for  all  e *.  and  for  all  n,  0 < n < N-l  . 
h h — 


We  will  prove  the  existence  and  uniqueness  of  the  solution  u^  of  Problem  (3.1);  for 


that  purpose  we  need  two  preliminary  lemmas. 
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For  each  n,  0 < n < N 1 , there  cuiitbponds  to  each  v t V a unique  function 

h h 

sfi  * y € #n  which  coincides  with  v on  Gn  . Let 

h h h h 

(3.2)  “ VV^V'  f°r  aU  Vh  ‘ Vh  ’ 

G 

£?n(vi  ,v  ) is  a bilinear  form  defined  on  V,  * V.  and  we  have 
n h h h 

Lemma  3 . 1 

For  all  n,  0 < n < N-l,  for  ail  u £ V,  and  all  v.e  V , 

- — h h h h 

(3.3)  8 n (u  , v ) = -((u  , -r: - v ) ) + ((grad  u , grad  v ))  + 

h h h dt  h n h h n 

G G 


, nil  n+1  , n n+0, 

+ (Uh  ' Vh  ‘n+1  ‘ lV  \ > n ' 


,,  ..  _n,  .11  . II  2 I,  n*li2  1 t n i2  1 i n+0  n .2 

(3.4)  8 (v  ,v  > = ||  grad  v ||  + j|  v , - j |v  | + j |v  - vh  | . 

G ft  ft  ft 

Proof:  Foi  each  v,  , the  function  - q^v.  is  continuous  on  G°  , hence 

h h h 

(•  ,tU)  = vnfi  and  formula  (3.3)  follows  at  once  from  (2.2).  Then,  we  deduce  (3.4)  by 
h h 

taking  v u , integrating  the  first  term  in  the  right  hand  side  of  (3.3)  with  respect 
h h 

to  t and  using  the  identity 

. n n*0.  1 | a | 2 1 i n+0 |2  1 i n+0  n ,2 

'v%  V ■ 2 I",  l0.  • ; K ln„  - 2 K - % In„  ■ 

We  will  also  need  the  following  elementary  lemma. 

Lemma  3. 2 

2 2 

Let  a,b,c  be  three  nonnegative  real  numbers  such  that  a b + ac  . Then, 
a <_  b + c . 

1 2 2 1/2  2 2 1/2 

Proof : We  have  a <_  — (c  + (c  + 4b  ) ) and  (c  + 4b  ) c + 2b  . ■ 

Theorem  3.1 


The  discrete  problem  (3.1)  admits  a unique  solution  u which  satisfies  the  estimate 
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2 +|lV2ni  l<lU°l  0 + .. 

G(o,t“)  2 h n"  n°  G(o,t") 


(3.5)  llgrad  ujj  _ n>  - 2,“h1„n  - 2l|“  >n0  " n 


)2  . 


for  all  n,  0 < n < N,  with  d = Max  {diam  il ( t ) } 

0 <_  t <_  T 

Proof:  For  each  n,  (3.1)  is  a system  of  linear  algebraic  equations  with  a square  matrix 

whose  order  is  equal  to  the  dimension  of  the  space  . Hence,  the  existence  of  the  solu- 

tion  u will  be  a consequence  of  its  uniqueness, 
h 

Taking  v>h  = in  usin9  <3-4><  summing  with  respect  to  n and  applying 

Poincard’s  inequality,  we  get 


2 . 1 i n , 2 

G(0,t“)  ~ “ ft” 


grad  V I..  ,n,  + 2 |Uh>  „n  - 


4lu°|2oH|f|l  n I!  uj  n< 

2 n°  G(0,t  ) n G(0,tn) 


iil“°|20  + dlifl!  n 

ft  G(0,t  ) 


grad  u. 


n"  G(0,tn) 


The  estimate  (3.5)  follows  at  once  by  application  of  Lemma  3.2  with 


a = ( ||  grad  u 


1 I n .2  . ) 
n.  + 2 |uh  ' n> 


1 |u°| 


G(0,t  ) 
and  c = d |1  f 


ft 


G(0,tn) 


The  uniqueness  of  is  a consequence  of  this  estimate. 

Remark  3.1:  The  equations  (3.1)  are  equivalent  to 

(3.6)  8 <uh'vh>  * ((f,vh)  £or  311  Vh  £ Vh  * 

T 

N-l 

with  (9(u  ,v  ) = 7 5n(u  ,v  ) . In  (3.6)  the  index  n does  not  appear.  However,  the 
h h „ h h 

n«0 

form  (3.1)  is  more  convenient. 

Remark  3.2.  For  each  n,  the  space  can  be  chosen  arbitrarily  with  no  relation  with 

the  spaces  for  v ? n;  in  particular,  the  spaces  need  not  have  the  same  dimen- 


sion for  different  values  of  n . 
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jt.. . 


1 

I 


Let  us  also  notice  that  unlike  other  methods,  the  present  method  does  not  require  a 

0 0 0 

preliminary  approximation  of  the  initial  function  u ; we  take  u^  «=  u 


4.  General  error  estimates 


In  this  section,  u is  the  solution  of  the  exact  problem  (2.1),  u^  is  the  solution 

of  the  discrete  problem  (3.1)  and  v,  is  an  arbitrary  function  in  the  finite  dimensional 

h 

space  V.  . We  estimate  the  error  u - u.  in  terms  of  u - v,  for  arbitrary  v c V.  . 
h h h h h 

In  later  applications,  we  will  choose  v,  equal  to  a certain  interpolating  function  of  u 

h 

(see  Section  5),-  thus  the  problem  of  estimating  u - is  reduced  to  the  problem  of 

estimating  the  interpolation  error. 

We  will  establish  the  following  result. 


Theorem  4 ■ 1 

Let  u be  the  solution  of  problem  (2.1)  and  u^  be  the  solution  of  the  discrete 
problem  (3.1) . Then, 


||  grad  (u-u  ) || 

G(0,t  ) 


/2 


|u"  - u"| 

h 


n-1 


v=0 


(4.1)  £ •d  d (£  ||  — (u-v^)  ||  2 v> 1/^2  + || grad  (u-v^) 

— G 


G(0,t  ) 


n-1 

_ I V Vi  r I V+0 

+2  Max  |u  - v | + 2 ) |v 

1 < v < n v=l 


for  all  functions  v,  e and  for  all  n,  1 < n < N 

h h 


For  the  proof  of  this  theorem  wc  will  need  the  following  lemma. 


Lemma  4 ■ 1 

Let  an  and  bn,  1 £ n N,  be  two  sequences  of  non-negative  real  numbers  which 
satisfy 

n-1 

(4.2)  (a")  + (bn)  <_  aan  + 6bn  + £ yVbV,  for  all  n , 

v— 1 

where  a, 6 and  yV  for  1 < v < N-1  are  non-negative  real  numbers.  Then 

(4.3)  a"  + bn  < /2  (a  + S + j yV)  . 

v=l 
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r 


Proof  of_Lemma  -1 . 1 
n 


Let  c 


n 2 n 2 1/2 
((a  ) t (b  ) ) 


Then,  (4.2)  yields 
"-1  v v 


(cV  l i a<-e)cn  + l yV  c°,  for  1 < n < N . 
v=l 


t 1 < n < N,  be  the  sequence  of  non-negative  real  numbers  which  satisfies 


(4.4) 


n-1 

(dn)2  *=  (a+g)d‘l  + l Y d’',  for  1 ± n £ N 
vel 


We  have  c1  < d1  = a + 8 and  d"  > n + 0 for  all  n . By  mathematical  induction  we 

prove  that  c"  < d"  for  all  n:  assume  cV  < d“  for  v=l,2 n-1;  then  g (c"l  < g (d") 

with  g (y)  = y2  - (ci+8)y;  since  g(y)  is  increasing  for  y > a + 5,  we  deduce  c < d . 
On  the  other  hand,  we  have  g(d")  <_g(dn+1)  for  all  n,  therefore  d"  < dn+1  for  all  n 


, v ,n 

and  (4.4)  yields  after  replacing  d by  d : 


. n-1 

n r V 

d <^a  + B+  2.  Y 
v=l 

Finally,  we  have  a"  t h"  < /2  c"  /I  d",  which  ends  the  proof  of  the  lemma. 
Proof  of  Theorem  4 . 1 

Taking  G = g"  and  <f  = in  (2.3)  and  subtracting  from  (3.1),  we  get 

S (U  - u_ ) =0,  V V>.  e t"  . 

' n h h n h 

G 


Hence,  we  have 

(4.5) 


S n(u-uh,u-q(r,)uh>  = ^n(u-uh,u-q(n)vh) 


for  all  v e V 
h h 


Let  tfn(u-u  , u-v  ) = & (u-uu,  u-qlv),  which  is  an  extension  of  the  definition 

h h ^.n  h h 

-» 

(3.2).  The  formulae  (3.3)  and  (3.4)  of  Gemma  3.1  are  also  valid  for  9 (u~lVu''Vh>  with 


un+0  . u11  « u(*,tn)  . Hence,  (4.5)  yields 
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i , 2 1 . n+1  n+1 1 2 li  n n .2 

||grad(u-uh)  11^  + j 'u  ' Uh  ^n+1  ' 2|u  'Uh  1 + 


1 | n+0  n 1 2 

-“J0„ 


= - ((u-u,  , — (u-v  )))  + ( (grad (u-u  ) , grad (u-v  )))  + 

h dt  h Gn  n n Gn 

, n-fl  n+1  n+1  n+1  . n n n n 

+ (u  - u , u - V ) , - (u  - u,  , u - v)  + 

“h  h nn+l  h h nn 

, n n n+0  n,  , _ , . „ . 

+ (u  - u,  v,  - v,  ) , for  0 < n < N-l  . 

h h h nn  — ~ 

By  summation  and  by  application  of  Poincare’s  inequality,  we  deduce: 
ii  ,,  , i|2  li  n ni2  1 nr^  | v+0  Vi2 


(4.6)  <_  II  grad  (u-u  ) ! 


f n-l 

n (d<n 

),t  ) ^ y=o 


,T<u-vJ  II2,,) 1/2  + llqradCu-v  ) 


at’  h'"Gv 


G (0 , t ) 


i n ni  I n ni  , r i v Vi  i v+0  Vi 

+ |u  -u|  |u  -v|  + l |u  -uh|  )vh  - vh)  , 

ft  ft  V=1  ft  M 


for  1 < n < N . 


Now,  wo  will  make  use  of  Lemma  4.1.  Lot 


a"  = || grad  (u-u  ) ||  , bn  = |Un  - u£|  , 

" G(0,  t ) J5 


a"  = d(  l ||  jTIu-v  ) ||  2V)1/2  + ||  grad  (u-v. 
v=0  G 


w ' 11 

h G(0,tn) 


_n  nr  I v vi  . n nr  i n+0  n 

B = /2  Max  |u  -v  | and  y = /2  |v  -v 


h h’nn 


Then,  (4.6)  yields 


, n.  2 n n „n,  n . r v 

(a)  a a + Bb  + I y b , for  1 £ n £ N . 
v=l 


Since  the  sequences  (an)  and  (b")  are  increasing,  we  deduce: 


(an)2  + <hn)2  < nlan  + bV’  + l y\^ , for  1 £ n < l < N 

v=l 
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Applying  l^emma  (4.1)  we  get: 


a"  + b"  < (a1  + Bl  + £ yV)  . for  1 £ n 1 1 


The  estimate  (4.1)  of  Theorem  4.1  follows  by  taking  n - J.  . 

Let  us  notice  that  formula  (4.5)  also  gives  an  estimate  for  the  expression 

"j1  |uV+°  - |2  which  involves  the  jumps  of  the  function  uh  . 

v=°  h h „v 

The  following  theorem  is  a variant  of  Theorem  4.1. 


Theorem  4 . 2 

Assume  that  Vh  n C°  ( i T>  is  non  empty  (where  C°  ( j.  T)  is  the  sgace  of  all  eontirv; 

nous  functions  defined  on  M ^ - Then,  we  have  the  estimate 


|| grad  (u-u  ) ||  + — 

h G(0,t  ) /2 


1 inn. 

+ — u -u. 

r-  ' h „ n 

✓2  n 


< /2  d II  ^ (u-vh) 


G(0,t  ) 


/2  ||grad(u-v 


hG(0.tn) 


_ i n n | 

+ 2'u  ' vh  jjh  ' 

for  all  functions  vh  c vh  n C° 1 * J and  for  all  n.  1 1 n < N . 

Proof:  In  this  case  we  have  |v^°  - v*|  = 0 . Then,  (4.6)  yields 

(a")2  + (t>n) 2 < a"an  + fiV,  where  a",  bn  and  a" 

are  defined  as  in  the  proof  of  Theorem  4.1  and  0 = ^2  |u"-v^n  . it  follows  antbn 

< /2(an  ♦ 6°)  . 

Notation : We  introduce  now  certain  general  notations  which  are  needed  in  the  following 


section. 


x x .....x  = space-coordinates, 

12  m 

j = (j  ,jx = multi-index  with  jQ,j1 3m  - ° 

(>l jm)  ' 

lil  = l ll,l*  IV I = l • 
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3ju  - a^u/at"*0  ax^1  ...  ax’*” 

L m 

aj’u  « ah' L/ax^1  ...  ax^ra  , 

° 'uf-s 

A 


DSu  ||  = ( I ||aju  ||^)  , for  any  integer  s > 0 and  any  G c J> 

*-  **  ' G T 


|DSu|n(t)  - < 


a-5  u(*,t)  lj}(t)>1/2'  for  any  0 £ 1 £ T 


The  two  foregoing  expressions  are  defined  provided 

u e H®CJ>  ) n C° ( [0 ,T]  ; HS(f5(t))) 


5 . Space - 1 i me  fjji  j t r < 1 ement*. 

In  this  section,  we  mak«  a specific  choice  for  the  spaces  we  use  space-time 

finite  elements  as  in  [ 3 ) , [4],  Ill],  112]  . For  simplicity,  we  assume  that  the  domain 
£ ^ is  polyhedral;  if  not,  we  should  use  curved  finite  elements  near  the  boundary, 

5.  a)  S i mp  1 i c i a 1 elements 

Let  h be  a "small”  positive  number  and  for  each  n,  0 < n < N-l,  let  J ” be  a 

— — h 

finite  set  of  (m+l)-simplices  K which  satisfy  the  conditions 

(5.1)  Gn  IlK;  Kf  l")  , 

h 

(5.2)  h (K)  <■  i,  for  all  K < 3 " , 

h 

where  h(K)  denotes  the  diameter  of  K , 
o o n o 

( ^ - 3 ) K n K'  if  , » K,K*  t Z . , where  K denotes  the 

h 

interior  of  K , 


(5.4) 


If  a vertex  of 
vei  tex  of  K'  , 


K belongs  to  K' , then  it  is  also  a 

for  all  K,  K ' e J . 

h 


Figure  2:  A triangulat ion  of  the  subdomains  cn  ^ and  in  the 

one-dimensional  case. 
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We  will  denote  7^  = { u 3^;  0 _<  n <_  N-l  } . Let  k be  an  integer,  k ^ 1;  let 

be  the  set  of  all  polynomials  of  degree  < k with  respect  to  the  variables  t,  x , ....  x 

~ In 


and  let  P (K)  be  the  set  of  their  restrictions  to  K . For  each  n,  we  choose  4>? 

K n 

equal  to  the  space  of  all  continuous  functions  defined  on  G n which  vanish  on  the  lateral 
boundary  of  Gn  and  whose  restriction  to  each  element  K coincides  with  a polynomial  of 
degree  < k,  i.e. 

(5-5)  *"  = {*.  ; # € $(g"),  ^ I £ P (K),  V K e 7 "}  . 

h h h h|K  k h 

For  the  statement  of  the  next  theorem,  we  need  the  following  definition.  Let 
denote  the  set  of  all  the  edges  of  a (m+l)-simplex  K,  with  m+1  = p > 2,  and  let 
8(D,D')  denote  the  angle  of  two  arbitrary  straight  lines  D and  D'  in  1R  p . 

Definition  5.1 

The  condition  angle  of  a p-simplex  K is  the  angle 


(5.6) 


6(K)  *=  Max  min  6(D,D*) 
De*P  D*eEK 


Remarks:  For  all  non- degenerate  p-simplices  K,  we  have  0 < 0 < x/2  . The  condition 

angle  0 ( K)  approaches  n/2  if  and  only  if  all  the  edges  of  K are  almost  parallel  to  a 
same  hyperplane.  In  the  case  p = 2,  K is  a triangle  and  6 ( K)  is  eqtial  to  one  half  of 
the  largest  angle  of  K . 

Theorem  S ■ 1 

Let  the  space  be  chosen  according  to  (5.5).  let  u be  the  solution  of  problem 

tot 

(2.1)  and  be  the  solution  of  problem  0 . 11  Assume  u e H ( j<  ) n C ( [0,T]  ; H (fl  (t)  ) ) 

where  t is  the  smallest  integer  such  that  t >_  k+1  and  t > (m+l)/2  . Assume  that  the 

discretization  of  J>  satisfies  the  conditions 
T 


(5.7) 

(5.8) 


h < 1,  Nh  < X 


0 (K)  <_  eQ  < »/2,  for  all  K t 7h  , 


where  X and  0Q  are  given  positive  constants.  Then, 


(5.9) 


||grad(n  -u)  ||  + |u"-un|  < y h 

G(0,t  ) n a 
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for  1 < n < N,  with 

t 

l 

s=kti 


= C }'  ( (cos  0Q) 


- 1 ||  s 

D u 


'ccct",  4XK<t»,Dttt''t,,0w’  ' 


where  C is  a positive  constant  which  depends  only  Oh  la,  k and  d . Moreover,  we  have 
the  following  estimate  for  the  jumps  of  u ^ : 

N-l 


(5.10) 


r | n+0  n|2,l/2  vk 

( l \%  ' V ) ' < C'Yh 

n = 0 


where  C‘  is  an  absolute  constant . 


Proof : For  each  element  K,  let  .0  (K)  be  the  set  of  the  points  of  K whose  barycentric 

coordinates  with  respect  to  the  vertices  of  K are  multiples  of  1/k  and  let  Tl^  de- 

0 - 

note  the  interpolation  operator  which  associates  to  each  function  f c C (K)  the  unique 
function  H c P ^ (K)  such  that  : v>  at  all  the  points  of  J (K)  (see  118]  , (5]). 

0 7T 

Let  II  denote  the  interpolation  operator  which  associates  to  each  function  c C w ) 
h T 

the  unique  function  II.  c V,  such  that  II  . = II  r y in  the  interior  of  K,  for  all 

* h h h K K 

K c Z.t  where  r denotes  the  operator  of  restriction  to  K . Note  that  the  interpolated 
h K 

function  nt  ^ is  in  general  discontinuous  at  the  time  t = tn  since  the  sets  J ^ 
h n 

and  IT  ™ are  independent. 

h 

We  will  use  the  estimate  of  Theorem  *1.1  with  v.  » n . u (note  that  II  ,u  is  defined 

n n n 

since  the  hypotheses  of  the  present  theorem  together  with  Sobolev’s  imbedding  theorem  im- 
ply th.  t u is  continuous  on  > ) . Let  Z * (resp.  «T  ' ) denote  the  set  of  all 

r 1 T h,  n h,n+0 

m-simpl ices  K'  which  lie  in  the  hyperplane  t - tn  and  which  are  a face  of  an  element 

Using  t inequality 


K c Z 1 (resp.  3 J1) 
n n 


, V+0  V,  i V+0  Vi 

|vh  * V - ,Vh  ~ U * 


I v vl 
|vh  - u I 


we  deduce  f roin  (4.1)  •- 


|grad(u-uh)||>T  + 7J  - 


(5.11) 


Kr  J , 


<^d(  l ||A  (u.nhU)||2,V2  + 


2 1/2 
K 


* ( l || grad ( u - nhu)||^)1^2  + 

Kt  J h 


-If.- 


+ 2 


N 1^<n  ( J . lU-nhulK,)]/2  + 

K 4 ^h.v 


+ 2 ( N- 1 ) Max  ( l |u  - n u|j\)1/2  , 

1<V<N-1  K’c  3 ' h K 

h,v+0 


v+0 

with  n u = (Thu)  on  each  K'  c 3 ' 

h h h , v+0 


(*) 

Theorem  2.3  of  [10]  yields: 


— (u-nhu)||K  + ||grad(u-nhu)||K 


(5.12) 

for  h < 1 and  for  all  K c 3 

(5.13) 


h * 


u - nhulK,  I c2h 


k+1 


C1  hk  f 

COS  0 l 

0 s-k+1 


l |nsu|K.  , 

s=k+l 


| DSU 


for  h < 1 and  for  all  K'  £ 3 u 3 ' , 

— h , v h , v-t  0 

where  and  are  constants  which  depend  only  on  k and  m . The  estimate  (5.9) 

follows  at  once  from  (5.11),  (5.12)  and  (5.13),  where  N can  be  replaced  by  an  arbitrary 
value  of  n . The  estimate  (5.10)  is  obtained  by  using  (5.9)  in  the  right  hand  side  mem- 
ber of  (4 . 6)  . ■ 

Remark  5.1:  The  condition  (5 . 8)  which  was  established  in  [ 10]  is  an  improvement  on  the  classical  con- 
dition of  Zlamal  [25]  and  Ciarlct-Raviart  [5].  It  is  valid  for  a general  class  of  interpo- 
lation operators  which  contains  the  operator  n„  considered  above  and  for  general  error 
estimates  in  the  Sobolev  spaces  p integer  0,  1 £ 9 i °°  • The  same  condition 

has  been  established  independently  from  the  author  by  BabuSka  and  Aziz  [1]  in  the  particular 
case  m = 1 (the  elements  K are  triangles),  k = 1 or  2,  W = 1 and  q = 2 . 

Remark  5.2:  Suppose  that  the  triangulation  3,  satisfies  the  condition 

h 


(5.14) 


h,n 


3 ’ . , for  1 < n < N-l  . 

h,n+0  — - 


0 ,~7 


Then,  nhu  e C (^  and  wc  can  apply  Theorem  4.2.  We  get  the  same  estimate  as  in  Theorem 


<*) 


There  is  a typographical  error  in  the  statement  of  this  theorem:  one  should  read 

k-m+l ,p  k + f,p  . , ,,k-m+l,p  . ,,k+l,p 

W and  W instead  of  W and  W . 
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5.1,  but  without  the  hypothesis  Nh  < X;  in  the  expression  of  y,  the  constant  X must 
be  replaced  by  1 (or  by  h ) . 

Remark  5.3:  The  condition  Nh  < X implies  that  the  average  value  of  the  time-steps 

should  not  be  too  small;  it  should  be  3t  least  of  the  same  order  as  h . In  particular, 
we  can  take  N - 1,  which  yields  a completely  implicit  method  for  the  whole  domain  J> 
The  advantage  of  dividing  the  interval  (0,Tj  into  sub-intervals  is  to  split  the  global 
system  of  discrete  equations  (1.1)  into  subsystems  each  of  which  corresponds  to  one  sub- 
domain  G . 


S . b ) Prismatic  elements 

For  simplicity,  we  will  use  the  vocabulary  corresponding  to  the  case  m = 2;  thus, 
we  will  say  "triangle"  instead  of  ”m  simplex",  "plane"  instead  of  "hyper-plane" , "prism" 
instead  of  "hyper-prism",  .. 


For  each  n,  0 < n < N-l,  let  3 . be  a finite  set  of  elements  K=K.,l<i<I, 
- h i — — 

which  satisfy  (5.1),  (5.2),  (5.3)  and  the  following  properties:  the  section  of  each  ele- 
ment by  the  plane  t = tn  (resp.  t = t"  ) is  a triangle  t"4^  (resp.  t"  + *)  and 

the  section  of  by  any  plane  inconstant  , tn  <_  t < tn  ^ , is  a triangle  whose  vertices 

are  located  on  three  straight  segments  [P^  p^4*]  where  p ^ 4 is  a vertex  of  Tn+^ 

and  P^  * is  a vertex  of  t'/4*;  the  element  Kn  is  a distorted  prism  (see  Figure  3); 
in  general,  it  is  not  a polyhedral;  the  triangles  Tn4t)  and  T*14*  are  called  the  bases 

of  the  element  Kn  and  their  vertices  are  called  the  vertices  of  yn  . All  the  vertices 
i i 

of  the  elements  K e 3 " must  satisfy  condition  (5.4).  The  sets  3 " and  3 n+^  are 
h h h 


independent;  in  particular,  the  set  3 1 of  all  triangles  Tri  + ^ which  are  a base  of 

h,n+0  i 

one  element  K c 3 !*  is  independent  of  the  set  3 of  all  trianqles  T1'  ^ which  are 

h h,n  l 

a base  of  one  element  K e 3 ^ (this  is  a difference  with  the  continuous  firite  element 

method  used  in  (41 ; in  that  method  we  imposed  3 ' = 3 ' ) (See  figure  4)  . 

h,n+0  h,n 

Let  k'  and  k"  be.  two  integers,  k'  > 1,  k"  > 1 . For  each  element  K « 3 " , 

— — n 

let  ( K)  be  the  space  of  all  functions  defined  in  K whose  restriction  to  each 

section  of  K by  a plane  t=constant  is  a polynomial  of  degree  k'  with  respect  to  the 
space  variables  and  whose  restriction  to  each  edge  [P^4°  P^*1)  is  a polynomial  of  degree 
k”  with  respect  to  t . We  take 
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(5.15) 


< - {V  *h  £ Vh  £ V,k"W'  K ‘ K] 


where  denotes  the  res'  iction  to  K . 


Let  k = min{k'  ,k"}  . t be  the  smallest  integer  such  that  t £ k + 1 and 

l > (m+l)/2  . Let  t‘  be  the  smallest  integer  such  that  1'  £ k'  +1  and  1'  > m/2  . 

For  each  K f l”,  let  h(K)  be  the  diameter  of  K,  h'  (K)  = (t*1  k - tn)  be  the  "height" 
h 

of  K and  p (K)  be  the  minimum  for  t"  < t < tn,k  of  the  diameter  of  the  largest  circle 

_ , _ n 

contained  in  the  triangular  section  of  K by  the  plane  t^constant.  Let  o ^ - tl  d ^ ; 

0 < n < N-l)  . Then,  the  following  result  holds: 


Theorem  5 ■ 2 

Let  the  space  be  chosen  according  to  (5.15).  Let  u be  the  solution  of  problem 

(2.1)  and  u.  bo  the  solution  of  problem  (3.1) . Assume  u ll?  (J^)  nC°([0,T);  H (h(t)))  . 
Assume 


(5.16)  Nh  £ X , 

(5.17)  li  (K)  £ h 1 1 , 

(5.18)  h'  (K) /h (K)  £ 0Q  , 

(5.19)  p(K)/h(K)  > Oj,  for  all  K f , 

where  A,  oq  and  are  positive  constants.  Then, 

(5.20)  ||grad(u  -u)  ||  + |u"  - u | n < Yh  + Y'h 

' G(0,t  ) ft 


for  1 < n < N,  with 


Y = C 


l l|DSu 


s=k+ 1 


G(0,tn) 


and 


Y1  *»  \C‘  l Max(|DSu(-,t)|ft(t);  0 £ t £ t")  . 


s-k'+l 


where 


C and  C'  are  two  positive  constants  which  depend  only  on  m,  k , k , 0^, 


d . Moreover 


(5.21) 


<Y  |u-°  - unJ2)1/2  £ C0(/  i Vhk) 

n-0 
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where  C is  an  absolute  constant. 

Q 

For  the  proof  of  this  theorem,  we  need  a preliminary  result  concerning  the  inter- 
polation error  in  each  of  the  elements  K . 

For  each  K e j”,  let  J(K]  bo  the  set  of  the  points  of  K which  are  located  in 

h 

one  of  the  planes  t = tn  + - tn)  with  j integer,  0 < j < k",  and  whose 

baryccntric  coordinates  with  respect  to  the  vertices  of  the  corresponding  triangular 
section  of  K are  multiples  of  1/k'  . Let  be  the  interpolation  operator  which 

associates  to  each  function  y e C^(K)  the  unique  function  Il^v'  £ Q^,  (K)  such  that 

II  * = ■p  at  all  the  points  of  J(K)  . Then,  the  following  result  holds. 

Lemma  5 . 1 

Suppose  that  K satisfies  conditions  (5.17),  (5.18)  and  (5.19).  Let  £ be  the 
smallest  integer  such  that  l >_  k + 1 and  £ > (m-i  l)/2  . Then,  for  all  functions 

t 

<t  € H (K)  , we  have 

. I 

(5.22)  ||d(*  - n *)||  <cxh  [ I|d‘v>|Ik  . 

s-k+1 

where  is  a constant  which  depends  on  m,  k',  k" , and  . 


Proof  of  Lemma  5.1 

Let  P (K)  be  the  space  of  the  restrictions  to  K of  all  polynomials  of  degree  £ k 

with  respect  to  the  variables  t,  x, , . ..,  x The  space  p (K)  is  invariant  under  the 
* 1 m k 

interpolation  operator  H . Hence,  applying  Lemma  2.5  of  [10]  which  is  a variant  of 
Lemma  7 of  Ciarlet-Raviart  “(5]  , we  get 

(5.23) 

where  C (K)  is  a constant  which  depends  on  the  element  K and  on  the  operator  H . 

1 K 

Let  us  assume  h(K)  = 1;  then,  we  have  P(K)  >_  0Q  and  h' (K)  ^ . The  set  \ of 

all  the  elements  K which  satisfy  these  conditions  is  compact.  Taking  the  maximum  of 
CjOO  for  all  K c V w can  replace  (K)  in  (5.23)  by  a constant  Cj  which  depends 
only  on  m,  k',  k“,  oQ  and  Oj  . The  estimate  (5.22)  for  h <_  1 follows  by  a simple 


|D(*  - I!K*>  )||K  < CjfK)  l I! II j.  , W £ n (K)  , 

s=k+l 


change  of  scale. 


Proof  of  Theorem  5 . ? 


Let  II  be  the  interpolation  operator  which  associates  to  each  function  v>  e ) 

n ip* 

the  function  II.  e V.  such  that  II  = II  r in  the  interior 
h h h K h 

of  K,  for  all  K c . Theorem  4.1  with  vh  = I^u  yields  (5.11)  as  in  the  proof  of 

Theorem  5.1.  The  estimate  (5.20)  follows  by  application  of  Lemma  5.1,  of  the  estimate 
(5.13)  with  k ' instead  of  k and  of  condition  (5.16).  Then,  the  estimate  (5.21)  follows 
from  (4.6). 

I 


6.  An  example 


In  this  section,  we  give  an  explicit  expression  of  the  discrete  equations  in  a 

simple  case:  uniform  rectangular  grid,  approximation  of  degree  1.  Thus  we  get  a finite 

difference  scheme  which  is  a particular  case  of  our  method.  We  will  assume  f = 0 . 

Let  Kn  , and  K?  be  two  neighbor  rectangular  elements  with  vertices  Pn  ,,  Pn+* , 
l-l  l l-l  l-l 

n n+1  n n+1  , n , n,  , n+1  n 

P.,  P.  , P..,,  P.^,  where  P.  = (x.,  t ),  x . , - x . = x . - x . , = Ax,  t - t = At  . 
l l i+l  i+l  l l l+l  l i l-l 

I«t  un  = u (Pn)  = uT(x.)  and  un+0  = u"+°(x.)  = lim  { u (x. ,tn+  e)  ; e > 0,  e ■+  0 } 
l hi  hi  i hi  hi 

In  each  rectangle  k”,  u^  is  linear  with  respect  to  each  variable  x and  t separately 
and  uniquely  determined  by  the  four  values  u^+^,  and  by  means  of  the  formula 


(6.1)  ii(x,t)  = (1  - x)(l  - t)  u + x(l  - t)  u ",  + (1  - x)t  u + x t u 


i+l  ' 


where  x = (x  - x/)/  Ax  and  t = (t  - t ) /At  . 


At  each  time  step,  the  values  u^  are  known  for  all  i and  we  must  determine  the 

values  u"+0  and  u^+^  • The  corresponding  equations  are  obtained  by  writing  the  integral 

relation  (3.1)  for  linearly  independent  test-functions  v>.  . For  example,  we  can  take,  for  each 

i,  the  two  test  functions  c*1'1'  and  such  that 

h h 

r i if  p = p" 

^l,1)  (P)  = < -1  if  P = p"+1 

0 at  all  other  grid-points, 


1 if  P = P. 


0 at  all  grid-points  P / P, 


Thus,  we  get  the  following  equations. 


,1,  n+0  n , 2,  n+0  n 1 n+0  n ,, 

<?(Ui+l  ' Ui+1>  + I(Ui  * V + 6(Ui-l  - Ui-1))/At  + 


. 1,,  n+1  _ n+1  n+1,  , n+0  , n+0  n+0,,  ,2 

+ 6 Ui+1  - 2 Ui  + Ui-1>  - (ui+l  ’ 2ui  + Vl))/(Ax)  = 0 

,1#  n+1  n+0  2,  n+1  n+0,  1,  n+1  n+0 

(6(ui+l  ' Ui+1>  + J(Ui  - Ui  > + 6(Vl  ' Vl),/ftt  " 

- (T<ui+l  - 2u?+0  + u^O)  + 2 (Un  + 1 . 2un+l  + un+l))/(Ax)2  _ Q 
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(7.1)  b)  Bu  = 0 on  E 

T 

0 

c)  u * u in  0(0)  , 

where  A is  a differential  operator  of  order  2y  with  respect  to  the  space  variables, 

V positive  integer,  and  B is  a boundary  operator. 

We  will  give  a variational  formulation  of  problem  (7.1)  which  is  a generalization  of 

the  integral  relation  (2.3).  For  each  t £ [0,T] , let  V (t)  be  a subspace  of  the  Sobolev 
li  2 

space  H (fi(t))  which  is  dense  in  L (f)(t) ) and  let  V ( t)  be  the  dual  space  of  V ( t ) . 

For  any  domain  G = GUq.t^)  = {(x,t);  x £ fl ( t ) , tq  < t < },  let 


E (G)  = L2(TQ,T  ; V(t))  , 

6 ' (G)  = V'(t))  = dual  space  of  6(G)  , 

8(G)  = {*;  <f  € e (G),  £ t'  (G)  } . 


denote  the  norm  in  £ (G)  defined  by 


IJ’liM 


3^  *||  q>1/2<  with  3 ' = (3X 3m)  . 


for  all  'fct  ( G)  . 


Let  ((*,*))  denote  the  duality  between  6(G)  and  6 ’ (G)  obtained  by  extension 
G 

of  the  inner  product  in  L2(G)  . 

We  have  6(G)  C l2(G)  c { ' (G)  and,  by  a lerata  of  Lions  [14],  8(G)  c ] ; 

L2<(Ht)))  . 
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Let  us  assume  Ac  £ ( £ (G)  , £ ' (G)  ) , i.e.  A is  a linear  continuous 

operator  from  £ (G)  into  £ ' C G)  and  let  us  define 

(7.2)  aG(i|',*)  = ( (A*  ,<f ) ) , for  all  i|i  end  <t  in  £ (G)  . 

a (•,•)  is  a bilinear  form  defined  and  continuous  on  £ (G)  * £ (G)  . 

U 

Assume  that  it  satisfies  the  two  following  properties. 


Uniform  continuity : 


|aG(*,*)  | < Ml|  *||  6 (G)  ||  v-  ||  e(G)  . 


for  all  * and  ^ in  £ (G)  and  for  all  G,  where  M is  a constant  indepen- 
dent of  G . 

Uniform  ooorcivity: 

(7.4)  aG<*,*)  + « ||  v"  ||G  > r>!I^H2e(G)  ' 

for  all  * c £ (G)  and  for  all  G,  where  <5  and  n are  two  constants  in- 
dependent of  G . In  fact,  we  will  consider  only  the  reduced  case  6=0 

B t 

obtained  by  a standard  change  of  variable  u ->  u e , g > 0,  in  (7.1)  . 


(7.5)  flGU,v)  = -(U,f|-))G  + aG(*,vO  + (*^)C(T  ) + 

_(,t"^)n(T0)  • 

8 (•,•)  is  a bilinear  form  defined  and  continuous  on  g(G)  * 3 (G) 

G 

Variational  problem:  Find  u c 8 (i<  T>  such  that  u ( • , 0)  = u®  and 

(7.6)  6 G{  u,*)  = ((f,?))G  , for  all  <e  c 8(G) 


and  for  all  G = GfTp,!^)  c j,  T . 

This  formulation  is  equivalent  to  the  formulation  of  Lions  [14]  . For 
the  existence  and  uniqueness  of  the  solution,  see  [14] . 

A standard  integration  by  parts  shows  that  (7.6)  is  a weak  form  of  the 
differential  equation  (7.1a).  As  for  the  boundary  conditions,  they  depend 
on  the  subspace  V(t)  . We  give  two  examples 


Example  7,1 

Let  Y (t)  = . Then,  the  solution  u of  problem  (7,6)  satisfies  the 

Dirichlet  boundary  conditions 

YjU  «=  0 on  , for  j = 0,1,..., y-1  , 

where  YjU(*»t)  = 3^u  ( • ,t)/3u''  denotes  the  trace  of  order  j of  the  function  u(‘,t)  on 
f(t)  . 


Example  7 . 2 


(X) 


Suppose  y = 1 and  A = A = Laplacian  operator.  For  each  t,  let  T (t)  be  a 
subset  of  r(t)  and  Y (t)  = W ; <fi  £ H1(«(t)),  = 0 on  T(1)(t))  . Let  l^l)  = 


•d) 


{ u r it ) ; 0 < t < T)  and  rT  = IT  - l T - 

satisfies  (in  a weak  sense)  the  boundary  condition: 


Then,  the  solution  u of  problem  (7.6) 


u “ 0 °n  I1^'  » 


, „ r (2) 

s u+t— =0  on  ) _ 
\j  3v  **  'P 


where  s^  denotes  the  outward  normal  speed  of  propagation  of  the  boundary  T(t)  and 


3u 

3v 


denotes  the  outward  normal  derivative  of  u . 


7 . b)  Discontinuous  approximations 

Our  method  of  approximation  is  based  on  the  fact  that  relation  (7.6)  makes  sense 

even  if  we  do  not  assume  ^ £ t (G)  since  ^ docs  not  appear  in  (7.6).  Therefore, 

ot  at 

this  relation  can  be  used  to  define  approximations  which  are  not  in  g (^  ) . 

Let  Gn  be  defined  as  in  Section  3 and  let  be,  for  each  n,  a finite 

h 

dimensional  subspace  of  g (Gn)  . Let  be  the  corresponding  space  defined  as  in 

Section  3;  is  ^ finite  dimensional  subspace  of  £ ( i , but  g ( Jr^);  if 

2 

v f V , the  function  t -*  v (»,t)  £ L (£l(t))  is  in  general  discontinuous  at  the  times 


t=t,  0 < n < N 


As  in  Section  3,  we  denote  v"  = v^(»,tn)  for  n = 0,1,..., N 

h h 


and  v^  ^ = lim{v.  (",tn  + e);  c > 0,  c -♦  0 ) for  n = 0,1,..., N-l  . The  discrete 

h h 

problem  is  formulated  exactly  as  in  Section  3. 


Discrete  problem:  Find  such  that  u^  = u°  and 
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f 




with  £D(v  ,v  ) = /p  (v  , q(n)v.)  . 

h h Gn  h h 

We  also  use  the  obvious  inequality  ||  ||  ^ <_  ||  v^|| 

Proof  of  Theorem  7.2 

Same  argument  as  for  Theorem  3.2. 

*££i  ications  to  finite  elements 

As  in  Section  5,  the  finite  dimensional  spaces 
space-time  finite  elements.  Then,  the  function  v 

h 

can  be  replaced  by  an  interpolate  of  the  solution  u 
estimated  by  using  standard  results  of  approximation 


£(G)  • 


can  be  defined  by  means  of 
in  the  right  hand  side  member  of  (7.9) 
and  the  interpolation  error  can  be 
theory. 


I 

l 
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